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^ ; Abstract 



In this paper we consider FRW cosmology in F(R, <j>) gravity. It is shown that in particular 
cases the bouncing behavior may appears in the model whereas the equation of state (EoS) pa- 
rameter may crosses the phantom divider. For the dynamical universe, quantitatively we also find 

43" 
Q H ' parameters in the model which satisfies two independent tests:the model independent Cosmological 



Redshift Drift (CRD) test and the type la supernova luminosity distances. 
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1. INTRODUCTION 



There are cosmological observations, such as Super- Nova la (SNIa), measurements of the 
cosmic microwave background (CMB) temperature fluctuations by the Wilkinson Microwave 
Anisotropy Probe (WMAP), the large scale red-shift data from the Sloan Digital Sky Survey 
(SDSS) and Chandra X-ray Observatory, that disclose some cross-checked information of our 
universe. Based on these observations the universe is very nearly spatially flat, and consists 
of approximately 75% dark energy (DE), 25% dust matter (cold dark matter plus baryons), 
and negligible radiation, in which the DE may drive the cosmic acceleration expansion l|-[4|]. 

One possible candidate for DE known as the cosmological constant which is given by 

;he vacuum energy with a constant equation of state may justify the late time acceleration 

5J-J9]. However, its value is so small ( of order 10~ 33 eV), in comparison with the Planck 

scale (10 19 GeV) that suffer from fine-tuning and the coincidence problems 10(. Instead, 

there are some other DE models are made by some exotic matter like phantom (field with 



negative energy) or some other scalar fields ll|, 112J. In general, there are many 



approaches 



to explain the origin of DE, which can be broadly classified into two classes 13|, ll4| . Either, 

through introduction of a more or less exotic form of matter such as phantom, quintessence 

15j or k-essence 16], or through modification of gravity such as brane- world models 17J, 



scalar tensor theories 



Gauss-Bonnet dark energy 181 ]. f(R) gravity [19|]-J27|, scalar-tensor theories 28 



generalized 



2jj, etc. 



Furthermore, according to the above observational data there exists the possibility that 
the EoS parameter has a dynamical behavior, evolving from larger than -1 (non-phantom 
phase) to less than -1 (phantom one ), namely, crosses -1 (the phantom divide) currently or 
in near future. In the framework of general relativity, a number of attempts to realize the 
phantom crossing have been made by many authors 30||- 32| . In the present work, we study 
the crossing behavior in the generalized scalar -tensor theory of gravity, F(R, </>). 

A bouncing universe which provides a possible alternative to the Big Bang singularity 
problem in standard cosmology has recently attracted a lot of interest in the field of string 



theory and modified gravity 33|, |34[ . In bouncing cosmology, within the framework of the 



standard FRW cosmology, the null energy condition (NEC) for a period of time around the 
bouncing point is violated . Moreover, after the bouncing when the universe enters into 
the hot Big Bang era, the EoS parameter in the universe crosses the phantom divide line 



35l . |36| . There are cosmological models in the framework of modified gravity which address 
both bouncing universe and crossing the barrier and again in here we discuss these in the 
framework of F(R, <fi) gravity. 

To understand the true nature of the driving force of the accelerating universe, mapping 
of the cosmic expansion is very crucial |37|. Here, we examine two observational tests in 
different redshift ranges to explain the expansion history of the universe 38|, |39| . The first 



probe we investigate is "Cosmological Redshift Drift" (CRD) test which maps the expansion 



of the universe 
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£ 



4l| (For a review see J42[-(45|) and measures the dynamics of the universe 



directly via the Hubble expansion factor. The second probe is the observations of the 
luminosity distance - redshift relation for SNIa that again verifies the late-time accelerated 
expansion of the universe. 

This work is arranged as follows. In the next section we present the F(R, <p) gravity. We 
drive the field equations and investigate analytically the conditions for the EoS parameter 
in the model to manifest phantom crossing and numerically analysis both bouncing and 
crossing the phantom line. Section three is devoted to the observational tests by including 
them in our analysis through driving CRD test and luminosity distance for the model. We 
also visit the Chevallier-Polarski-Linder (CPL) parametrization model 46j in a comparison 
with our model and observational data. We finally conclude with a summary and remarks 
in section four. 



2. THE MODEL 



We start with the scalar-tensor theory of gravity where the scalar curvature R in the 
action is replaced by F(R, </>), 



S= I d A Xy/^jF(R,, 



(1) 



with g is the determinant of the metric tensor g^ v and is a scalar field. It may originally 
turns out that F(R, <fi) model which may be considered as a generalisation of F(R) gravity, 
(with or without matter) can be used to solve some of the problems in cosmology. However, 
in general, the construction of F(R, <fi) is not explicit and it is necessary to solve the second 
order differential and algebraic equations. In here, we consider a polynomial form of F(R, 0) 
as Yl^=o Pn{4>)R n , as an extension of the formulation given in 27| 47|, where the scalar field 



is without kinetic term. 

The action Eq. ([I]) is then rewritten as 

/oo 
d 4 Xv ^(J2Pn^)R n )- (2) 

n=0 

and its variation with respect to the metric tensor g^ gives us the field equations as 

-.00 00 

= -2^^ P «^) jRn + ^ + ^ n -V M V ! ,))^nP n (0)/?^ 1 . (3) 

n=0 n=0 

Moreover, variation of the action ([2]) with respect to the scalar field <fi gives 

00 

J2K(m n = Q, (4) 

n=0 

where i^(0) = ji . The field equations ([3]) corresponding to standard spatially-flat FRW 

universe for the time and spacial components of the metric yields, 

00 00 ,00 

= -?>H 2 ^nP n {(l>)R n - 1 -3HY,nP n {<P)R n ^ + 3H—^nP n (<p)R n ~ l 

n=0 n=0 n=0 



_. 00 
+ -£P B (0iT, (5) 

n=0 

= -( 3 H 2 + 2H)J2nRMR n - 1 + (^ + 2H- + HJ^nPni^R"- 1 

n=0 ^ ' n=0 

1 00 

+ ^Y, R M Rn - ( 6 ) 

n=0 

In comparison with the standard Friedman equations, ^pf = 3H 2 and ^^ = —3H 2 — 2H one 
can consider the model as standard model with the effect of the F(R, (f>) gravity modification 
is contributed in the energy density and pressure of the Friedman equations. After some 
algebraic calculation, one can read the equivalently effective energy density and pressure 
from the above equations as, 

— M 2 ( rl °° -,00 "| 

^n=0 ™WV ^ n=0 n =0 J 

7M"2 f / J2 J \ °° -,00 ^ 

p «" = a^V' 1 1* 5 + 2/f * + *) § nPMR "~ l + 2 g p " (0)fl ") (8) 

Using Eqs. ((Tj) and ([5}, the conservation equation can be read as, 

Peff + 3#p e //(l + w e// ) = 0, (9) 



where u; e // = ^^-. From Eq. (jSJ) one also obtain, 



H 




H 



Sjn B 



(10) 



where 



n=0 
oo 

B = ln^P n (0)iT. 



(11) 

(12) 



n=0 



To study the behavior of the scale factor of the universe and Hubble parameter, for simplicity, 
we consider two linear and quadratic forms of the polynomial F(R, 0). Numerically, we solve 
the equations (jlj)-(jHJ) for these cases and obtain the results given in Figs. (1) and (2). As 
can be seen, in Fig. (1), in linear case, for exponential behavior of the Pq{4>) and power law 
behavior of the Pi{4>), the scale factor bounces at t = 0. In Fig. (2), with similar exponential 
behavior of the Pq(4>) but in quadratic form of F(R, 4>) where P2((f>) is power law and Pi{4>) 
is logarithmic, we obtain different dynamics for the scale factor. In this case the universe 
has an inflection at t — 0. 
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Fig. 1: The graph of a and H plotted as functions of time for linear 

F(R,<f>) = Pi((j))R + P (<t>), Pi(0) = a , P (<t>) = e w , for a = -0.125 and = 

Initial value is 0o(O) = 0.5. 
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Fig. 2: The graph of a and H plotted as functions of time for quadratic 

F{R, 0) = P 2 (<P)R 2 + PMR + P o (0), A(0) = hi(0), Pi(0) = r and P o (0) = e^, 

for a = —0.125 and f3 = —0.5. Initial value is 0(0) = 0.5. 

At this stage we investigate the cosmological evolution of EoS parameter, and show 
that there are conditions that cause the EoS parameter crosses the phantom divide line. 
Analytically, one requires H = —\{p e ff + Peff) vanishes for a particular <fr and changes 
the sign after the crossing. One may explore this possibility by checking the condition 



ji 



(Peff +Pe//) 7^ when u> e ff —7- —1. From Eqs. ([5]) and (EJ), one gets, 



H = ~\\ 



(13) 



Also we have ^(p e // + Peff) = —2H or, 



H = -\{a+(2A- H)A - HA} 



(14) 



One can find that in order to have phantom crossing one of the following conditions might 
be satisfied when u e ff — > — 1: 

• (a) If A = or 2A = H then A ^ 0, 

• (b) If A = then A ^ and 2A ^ H. 

Numerically, As shown in Fig. (3), for both linear and quadratic F(R, 0), the EoS parameter 
crosses the boundary at t > in the past. 
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Fig. 3: The graph of w e //, plotted as function of time for linear and quadratic forms: 
Left panel: F(R, 0) = Pi(0)P + P o (0), Right panel: F(R, 0) = P 2 (0)P 2 + P\{<f>)R + P o (0), 

P 2 (0) = ln(0), Pi(0) = a and P o (0) = e^. In both, a = -0.125, (5 = -0.5 and 0(0) = 0.5. 

Similar to the difficulties arise in numerically solving the forth order nonlinear differential 
equations in F(R) gravity for the scale factor a(t) 48J, here in F(R, 0), even in relatively 
simple cases of linear, quadratic and cubic forms of the model, a large uncertainties on the 
boundary conditions have to be set to find out the scale factor. In a different approach 
in F(R) gravity, the authors in 48|], instead of choosing a parameterized expression for 
F(R) and numerically solving the equations with the given boundary conditions, relate the 
present day values of F(R) derivatives to the cosmographic parameters (i.e. the deceleration 
go, the jerk j , the snap Sq and the lerk l parameters) in order to constrain them in a model 
independent way and suggest the kind of F(R) theory could be able to fit the observed 
Hubble diagram. A similar approach which is beyond the scope of this work can be applied 
in our F(R, 0) gravity model, in a separate study. 



3. COSMOLOGICAL TESTS 



3.1. Cosmological Redshift-Drift Test(CRD) 



In this section, to detect the cosmic acceleration more directly, following 40J], 49|- 52|. 
we employ the time derivative of the redshift in terms of the Hubble parameter as 



z = H [l + z 



H n 



(15) 



where z measures the rate of expansion of the universe: z > (< 0) indicates the accelerated 
(decelerated) expansion of the universe, respectively. The redshift variation is related to the 
apparent velocity drift of the source: 

(16) 



l + z' 



where v = -J^- and H = lOOhKm/sec/Mpc [53|. Since the detection of the signals 
correspond to the redshift change and consequently velocity drift is yery weak, to measure 
the signals, observation of the LYa forest in the QSO spectrum 42]- 45| for a decade 
might be needed . In near future this requirement can be provided by a new generation of 
Extremely Large Tele scope ( ELT) which are equipped with aspectrograph to measure such 



a small cosmic signal 41] 



50 



Ml- 



In the following, we use three sets of data (8 points) for redshift drift generated by 
performed Monte Carlo. We first revise the CPL parametrization model and compare the 
F(R,(p) model with CPL model and observational data 40[], 50j-[52|. 

CPL model 

One popular parametrization, which explains evolution of dark energy is the CPL model in 
which in a flat universe the time varying EoS parameter in term of redshift z is parameterized 
by, 



uj(z) =w + wi(— — ). 
l + z 

In the model, the Hubble parameter also is given by 

H(z) 2 



H 2 



n m (i + zY + (i 



i L rn _ ) 6 



Sfo^dz 



where using Eq. (1171) we can obtain the following equation for Hubble parameter 



H(z 



HI 



}- = Q m (i + zf + (i - n ro )(i + ^a+-o+-i) x exp 



-3wi(tt-) 



(17) 



(18) 



(19) 



l + z- 

ues of u and uj\. The velocity 



In CPL model the parametrization is fitted for different va 

drift with respect to the source redshift is shown in Fig. (5) 40|. A comparison of the model 

with the observational data shows that the best fit values are for uq = —2.2 and uj\ = 3.5. 

F(R, 0) Model 

From numerical calculation, in Fig. (4), it can be seen that for three linear, quadratic 
and cubic form of F(R,4>), the EoS parameter crosses —1 at different z. In linear case it 
crosses twice in the past while in quadratic and cubic cases it crosses only once. It shows 



that for cubic form, the phantom crossing is more appropriate with the observational data 
[? ] as sometimes in the past the EoS parameter becomes zero and positive while this is 
not the case in linear and quadratic forms. In Cubic case, the EoS parameter crosses the 
boundary at about z = 1.4 and its present value is about —1.3 




Fig. 4: The graphs of w, plotted as function of redshift z. 

I: F(R, 0) = Pi(0)P + P o (0), 0(0) = 0.05, II: F(R, 0) = P 2 (0)P 2 + Pi(0)P + P o (0), 0(0) = 0.5 

III: F(R, 0) = P 3 (0)P 3 + P 2 (0)P 2 + Pi(0)P + P o (0), 0(0) = 0.5. P 3 (0) = sin(0), P 2 (0) = ln(0), 

Pi(0) = a and Po(0) = e^\ In all cases we assume a = —0.09 and /3 = —0.6. 

From equations (jSJ), ([6]) and using ,; - = — (1 + z)H(z) , z - , we obtain the following 



expression for EoS parameter, 



w 



e// 



(! + *)£ 
3r 



(20) 



where r = |L. Solving equation ( 1201) for H (z) we obtain, 

■H(z)i2 r n 



^0 



exp 



w, 



e//^ 



1 + 5 



-d5 



(21) 



where from numerical computation in our model we find w e //, and then by substitution into 
equation (|2T|) one obtains H(z) for our model. Then, using equation (Tl6|) . one calculates the 
velocity drift. Figure (5) shows the velocity drift against z for linear, quadratic and cubic 
F(R, 0). The figure also display the parameterized CPL model result 54j. In comparison 
with the CPL model and experimental data, our model is in a better agreement for cubic 
form of F(R, 0). 



0.5 



X) 




Fig. 5: The graph of v, plotted as function of redshift z. 

I: F(R, 0) = Px(0)P + Poit), 0(0) = 0.05, II: F{R, 0) = P 2 (0)P 2 + P x (0)P + P o (0), 0(0) = 0.5 

III: F(R, 0) = P 3 (0)P 3 + P 2 (0)P 2 + P(0)P + Po(0), 0(0) = 0.5. P 3 (0) = sin(0), P 2 (0) = ln(0), 

Pi(0) = a and Po(0) = e^^. In all cases we assume a = —0.09 and j3 = —0.6. 



3.2. Luminosity di,(z) and the difference in the distance modulus /i(z) 

Luminosity distance quantity, di{z), determines dark energy density from observations. 
In CPL model with H(z) given by (TT9|) . it is in the form of, 

1 + z f* dz' 



driz) 



H n 



\/n m (i + z 1 ) 3 + (i - n m )(i + 2 ')3(i+wo+"i) x J-^4^\ 



(22) 



where Q m is the matter density. In addition, the difference fi(z) between the absolute and 
apparent luminosity of a distance object is given by, fj,(z) = 25 + 5log w dL(z) . In Fig. 6, 
we compare our model with CPL parametrization model for luminosity distance and with 
the observational data for fJ,(z). It shows that the cubic form of our model is in a better 
agreement with the CPL model for luminosity. It also shows that again cubic form of the 
model fits the observational data. 
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n =0.25 





Fig. 6: The graph of luminosity di and the difference in the distance modulus fi, 

plotted as function of z for CPL parametrization and observational data for our model. 

I: F(R, 0) = Pi(0)P + P o (0), 0(0) = 0.05, II: F(R, 0) = P 2 (0)P 2 + Pi(0)P + P o (0), 0(0) = 0.5 

III: F(R, 0) = P 3 (0)P 3 + P 2 (0)P 2 + P(0)P + Po(0), 0(0) = 0.5. P 3 (0) = sin(0), P 2 (0) = ln(0), 

Pi(0) = a and Po(0) = e 13 ^. In all cases we assume a = —0.09 and (3 = —0.6. 

4. CONCLUSION 

In this paper, we study the evolution of the gravitational and scalar fields in modified 
gravity in which a reconstructed F(R, 0) as a polynomial function of R of degree n is 
introduced. We find that the behavior of the scale factor of the universe depends on the 
degree of the polynomial. For linear F(R, 0) a bouncing behavior and for a quadratic one 
an inflection behavior for the universe is shown. The evolution of the cosmological EoS 
parameter against time in both linear and quadratic cases, shows a crossing of the phantom 
divide line at t > in the past with a transition start from u e ff < — 1 to the u e ff > — 1 in 
linear case and from u e ff > — 1 to the u e ff < —1 in quadratic case. A comparison of the 
phantom crossing in terms of redshift z for linear, quadratic and cubic forms are given. 

With respect to the observational data it shows that the quadratic and cubic forms are 
in better agreement with the data. It shows that in both cases the current value of the 
EoS parameter is between —1 and —1.4, the crossing occurs at about z = 0.5 for quadratic 
and about z — 1.4 for cubic case, and transition from matter to dark energy dominated era 
occurs only in cubic case 
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We then analyze the model with respect to the CPL model and observational data. The 
first cosmological test is the CRD test. The variation of velocity drift for linear, quadratic 
and cubic forms of the polynomial with redshift z is shown in Fig. (5). Checking between 
our model with CPL parametrization model and observational data shows that cubic form 
is in better agreement with the experimental data. The second cosmological test is luminos- 
ity distance d^ and the difference in the distance modulus /i that the former is calculated 
in comparison with the CPL parametrization model and the later with respect to the ob- 
servational data. Again, for luminosity distance it shows that the the cubic from of our 
model is in better agreement with the CPL model. For fi(z), it show that again the cubic 
form of the model very well fit with the observational data. Although, due to numerically 
difficult computation we have not considered the higher order polynomial forms of F(R, 4>), 
but from these three linear, quadratic and cubic forms we may conclude that the presence 
of the higher order terms in the model better fits the model with the observational data. 
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